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Abstract 

Convex relaxations of the power flow equations and, in particular, the Semi-Definite Programming 
(SDP), Second-Order Cone (SOC), and Convex DistFlow (CDF) relaxations, have attracted significant 
interest in recent years. Thus far, studies of the CDF model and its connection to the other relaxations 
have been limited to power distribution systems, which omit several parameters necessary for modeling 
transmission systems. To increase the applicability of the CDF relaxation, this paper develops an ex¬ 
tended CDF model that is suitable for transmission systems by incorporating bus shunts, line charging, 
and transformers. Additionally, a theoretical result shows that the established equivalence of the SOC 
and CDF models for distribution systems also holds in this transmission system extension. 


Nomenclature 

N - The set of nodes in the network 
E - The set of from edges in the network 
E r - The set of to edges in the network 
i - imaginary number constant 
I - AC current 
S = p + iq - AC power 

V = v/-6 - AC voltage 

Z = r + ix - Line impedance 

Y = g + ib - Line admittance 

T = tZ6* - Transformer properties 
Y s = g s + ib s - Bus shunt admittance 
W - Product of two AC voltages 
L - Current magnitude squared, |/| 2 
6 C - Line charging 


s” - Line apparent power thermal limit 

9 a - Phase angle difference limit 

S d — p d + iq d - AC power demand 

S 9 — p 9 + iq 9 - AC power generation 

co,ci,C2 - Generation cost coefficients 

K(-) - Real component of a complex number 

$>(•) - Imaginary component of a complex number 

(■)* - Conjugate of a complex number 

| ■ | - Magnitude of a complex number, Z 2 -norm 

Z - Angle of a complex number 

x u - Upper bound of x 

x l - Lower bound of x 

x - Convex envelope of x 

x - A constant value 
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1 Introduction 


Convex relaxations of the power flow equations have attracted significant interest in recent years. They 
include the Semi-Definite Programming (SDP) [1], Second-Order Cone (SOC) [12], Convex-DistFlow 
(CDF) [6], and the recent Quadratic Convex (QC) [10] and Moment-Based [19, 20] relaxations. Much of 
the excitement underlying this line of research comes from the fact that the SDP relaxation has shown 
to be tight on a variety of case studies [14], opening a new avenue for accurate, reliable, and efficient 
solutions to a variety of power system applications. Indeed, industrial-strength optimization tools (e.g., 
Gurobi [9], cplex [11], Mosek [24]) are now available to solve various classes of convex optimization 
problems. In the context of power distribution systems, the relationships between the SDP, SOC, and 
CDF relaxations is now well-understood [16, 17]. In particular, the SOC and CDF relaxations are known 
to be equivalent [23, 3] and the QC and SDP relaxations are at least as strong as both of these [22, 5]. 

Thus far, studies of the CDF model and its connection to the other relaxations have been limited to 
power distribution systems, which omit several parameters necessary for modeling transmission systems 
(i.e. bus shunts, line charging, and transformers). This paper develops an extended CDF model that 
incorporates all of the parameters necessary to use the CDF model on transmission system test cases. 
The main contributions of this work can be summarized as: 

1. Developing a novel methodology for derivation of the classic CDF model. 

2. Utilizing that methodology to produce an extended CDF model that is applicable to transmission 
systems. 

3. Proving that the equivalence of the SOC and CDF models carries over into the transmission system 
case. 

The rest of the paper is organized as follows. Section 2 reviews the formulation of the simplest AC 
power flow feasibility problem (AC-PF). Section 3 develops a number of general properties of AC-PF 
that are utilized throughout the paper. Section 4 reviews the two well known relaxations SOC and CDF 
for the AC-PF problem. Section 5 is the heart of the paper and extends all of the previous results to 
an extended AC-PF formulation incorporating bus shunts, line charging, and transformers. Lastly, the 
paper concludes with Section 6. 


2 AC Power Flow 

This section reviews the specification of the AC Power Flow (AC-PF) feasibility problem and introduces 
the notations used in the paper. In the equations, bold face is used to reflect constants while script 
indicates variables. Capital letters refer to complex numbers while lower case are real numbers. 

A power network is composed of a variety of components such as buses, lines, generators, and loads. 
The network can be interpreted as a graph ( N , E) where the set of buses N represent the nodes and 
the set of lines E represent the edges. Note that I? is a set of directed arcs and E R will be used to 
indicate those arcs in the reverse direction. To break numerical symmetries in the model and to allow 
easy comparison of solutions, a reference node r e N is also specified. 

The AC power flow equations are based on complex quantities for current I, voltage V, admittance 
Y, and power S, which are linked by the physical properties of Kirchhoff’s Current Law (KCL), i.e., 

i?-it= H hi (i) 

( i,j)eEuE R 

Ohm’s Law, i.e., 

hi = Yij(Vi — Vj) (2) 

and the definition of AC power, i.e., 

Sij = VJ* (3) 
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Combining these three properties yields the AC Power Flow equations, i.e., 

S?-S? = £ yieN (4a) 

(.i,j)eEvE R 

Sii = Y* | V| 2 - Y*jViV* (t, j)eE u E R (4b) 

These non-convex nonlinear equations define how power flows in the network and are a core building block 
in many power system applications. However, practical applications typically include various operational 
side constraints on the flow of power. We now review some of the most significant ones. 

Generator Capabilities AC generators have limitations on the amount of active and reactive power 
they can produce S 9 , which is characterized by a generation capability curve [13]. Such curves typically 
define nonlinear convex regions which are often approximated by boxes in AC transmission system test 
cases, i.e., 

Sf sS S? Sf u VicN (5) 

Line Thermal Limit AC power lines have thermal limits [13] to prevent lines from sagging and 
automatic protection devices from activating. These limits are typically given in Volt Amp units and 
constrain the apparent power flows on the lines, i.e., 

IS.,10" V(i,i)£EuE R (6) 

Bus Voltage Limits Voltages in AC power systems should not vary too far (typically +10%) from 
some nominal base value [13]. This is accomplished by putting bounds on the voltage magnitudes, i.e., 

v\ < |V| Vie IV (7) 

Phase Angle Differences Small phase angle differences are also a design imperative in AC power 
systems [13] and it has been suggested that phase angle differences are typically less than 10 degrees in 
practice [21]. These constraints have not typically been incorporated in AC transmission test cases [26]. 
However, recent work [4, 10] have observed that incorporating Phase Angle Difference (PAD) constraints, 
i.e., 

- 0$ sc Z(VV/) ^ 0$ V(i,j) e E (8) 

is useful in the convexification of the AC power flow equations. For simplicity, this paper assumes that 
the phase angle difference bounds are symmetrical and within the range (—7r/2, 7 t/ 2), i.e., 

0 Z 0# Z | (i,j)eE 

but the results presented here can be extended to more general cases. Observe also that the PAD 
constraints (8) can be implemented as a linear relation of the real and imaginary components of ViV* 
[18], i.e., 

tan(-0£)» {ViV*) sg {W*) ^ tan(0^)3? {W*) V(i,j) e E (9) 

and that equation (4b) can be used to express these in terms of the S variables as follows, 

ViV* = I Vi I 2 - Z+Sii (10) 

These equations combined with (9) implement the PAD constraints in terms of the V and S variables as 
follows, 

tan(—0£)H (|V| 2 - ZtjSij) V 9 (|V| 2 - Z*^) + tan(o£)» (|V| 2 - Z*S tj ) V(t, j) e E (11) 
The usefulness of these alternate formulations will be apparent later in the paper. 
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Model 1 The AC Power Flow Feasibility Problem (AC-PF) 


variables: 


Sf VieN 
Vi Vi e N 
Sij (i,j)eEvE R 


subject to: 


v\ < \Vt\ =% vf Vi e N 


sf V Sf V Sf u VieN 
I V(i,j) e E v E R 
Sf -Sf= Y S V e N 

U „'\c 177, , J? R 


Z V r = 0 


(12a) 

(12b) 

(12c) 

(12d) 

(12e) 


( i,j)eEuE 



( 121 ) 

(12g) 


Other Constraints Other line flow constraints have been proposed, such as, active power limits 
and voltage difference limits [14, 18]. However, we do not consider them here since, to the best of our 
knowledge, test cases incorporating these constraints are not readily available. 

The AC Power Flow Feasibility Problem Combining all of these constraints yields the AC 
Power Flow Feasibility presented in Model 1. Constraint (12a) sets the reference angle, to eliminate 
numerical symmetries. Constraints (12e) capture KCL and constraints (12f) capture Ohm’s Law. Con¬ 
straints (12b) and (12g) capture the voltage and phase angle difference operational constraints. Finally 
constraints (12c) and (12d) enforce the generator output and line flow limits. Notice that this is a non- 
convex nonlinear satisfaction problem due to the product of voltage variables, ViV*, and is NP-Hard in 
general [25, 15]. 

3 Generic Properties of AC Power Flows 

Before deriving the various relaxations of the AC power flow equations, this section develops a collection 
of general properties of the AC power flows (AC-PF), which are utilized throughout the rest of the paper. 

Absolute Square of Ohm’s Law From definition (2), the absolute square of Ohm’s law is, 



(13a) 

(13b) 


Absolute Square of AC Power From definition (3), the absolute square of the definition of ac 
power is, 


s,jS* = ( Vii*yy?iii) 
l ^| 2 = |b»| 2 |fy | 2 


(14a) 

(14b) 
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Absolute Square of Voltage Products The absolute square of the voltage product is, 


(ViV*)(ViV*)* = (W*)(V*Vj) (15a) 

\ Vi V?\ 2 = \Vif\Vi\ 2 (15b) 

Line Loss From definition (4b), the power loss on line (i,j) is, 

Su + S jt = Y*(\Vi\ a - ViV* - V*Vj + \Vj\ 2 ) (16a) 

Voltage Drop From definition (4b), the so-called voltage drop property is derived by subtracting the 
power on each end of a line (i,j), solving for the \Vi\ 2 , \Vj\ 2 variables, and eliminating the Sji term. 

Sij - S jt = Y*m 2 - VV* + V*Vj - \Vj\ 2 ) (17a) 

|V| 2 - \Vj\ 2 = ZfjSii - Z*jSji + ViV* - V*Vj (17b) 

\Vi\ 2 - |V,| 2 = Z+Sij - Z*(Y* (|V| 2 - ViV* - V*Vj + |V,| 2 ) - Sij) + ViV* - V*Vj (17c) 

|V| 2 - |V)| 2 = ZfjSij - (|V| 2 - ViV* - V*Vj + |^| 2 ) + ZfjSij + ViV* - V*Vj (17d) 

|V| 2 - |F,| 2 = ZfjSij - (|V| 2 - ViV* - V*Vj + \Vj\ 2 ) + ( \Vi \ 2 - ViV*) + W* - V*Vj (17e) 

|V| 2 - l^l 2 = (ZfjSij + ZijSfj) - (|V| 2 - ViV* - V*Vj + IV,I 2 ) (17f) 


Equivalence of Line Flow Formulations A key observation from these general power flow 
properties is that the line loss and voltage drop properties provide an alternate formulation of the power 
flow constraints (4b), namely, 

(4b) ^ (16a), (17f) 

which follows from the fact that these two sets of constraints are simply linear combinations of one 
another. 


4 Derivation of the Relaxations 

This section derives two well known convex relaxations of Model 1, the SOC relaxation [12] and the 
Convex DistFlow relaxation [6] and reviews their equivalence. Although none of the results in this section 
are new, the derivation of the models presents a novel and systematic methodology that is leveraged in 
subsequent sections for developing the extended Convex DistFlow relaxation for transmission systems. 

4.1 The SOC Relaxation 

The SOC relaxation was first proposed in [12] and utilizes two key insights. First, by lifting the product 
of voltage variables ViV* into a higher dimensional space (i.e. the IV-space), 

Wi = \V\ 2 i e N (18a) 

Wij =ViV* V(t,j)e£ (18b) 

a linear relaxation of (AC-PF) is obtained. Second, the absolute square of voltage products property 
(15b) can be used to strengthen this IV-space relaxation, as follows, 

\Wij\ 2 = WiWj V(i,j) e E (19a) 

\Wij\ 2 H WiWj V(i,j) eE (19b) 

Notice that constraint (19b) is a convex second-order cone constraint, which is widely supported by 
industrial strength convex optimization tools (e.g., Gurobi [9], CPlex [11], Mosek [24]). 

The complete SOC relaxation of (AC-PF) is presented in Model 2 (SOC-PF). The constraints for 
KCL (12e), generator output limits (12c), and line flow limits (12d) are identical to the (AC-PF) model. 
Constraints (20c)-(20d) capture line power flow in the W-space. Constraints (20a) and (20b) capture 
the voltage and phase angle difference operational constraints. Finally constraints (20e) strengthen the 
relaxation with voltage product second-order cone constraint. 
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Model 2 The SOC Relaxation of AC Power Flow (SOC-PF) 


variables: 

S? Vie N 
Wt VieN 
Wij V(i,j) e E 
Sij V(i,j) e£u E r 

subject to: 

(12c), (12d), (12e) 


{v l t ) 2 Wi < «) 2 Vie N (20a) 

tan(-0£)ft (W«) 3 (Wy) tan(0£)ft (Wy) V(i,j) e E (20b) 

Sij = Y*W. i - Y*\l,j (i, j) e E (20c) 

Sji = Y*Wj - Y*W* (*, j) e E (20d) 

|Wy| 2 < WiWj V(i,j)eE (20e) 


4.2 The Convex DistFlow Relaxation 

The DistFlow (DF) model is a non-convex power flow model originally developed in [2] and later the 
Convex DistFlow (CDF) model was proposed in [6]. Both DF and CDF were originally designed for 
radial topologies, however, [7, 8] show that they can be interpreted as a phase-angle relaxation of the 
meshed AC power flow and hence are applicable to the study of meshed power networks. This section 
presents a novel derivation of the DF model and its relaxation from first principles. The value of this 
derivation is to establish a clear connection between (AC-PF) and the CDF model. 

The DistFlow model [2] can be derived utilizing three key insights. First, replace the line flow 

equations (4b) with their alternate formulation based on line loss (16a) and voltage drop (17f). Second, 

lift the model into the space of bus voltages and line currents (i.e. the L-space), 

Wi = \Vi \ 2 is N (21a) 

Lij = |^| 2 (|F i | 2 - VW* - V*Vj + |V)| 2 ) V(i,j) e E (21b) 

Third, use the absolute square of power property (14b) to strengthen the L-space relaxation, as follows, 

\S iJ \ 2 = WiL ij V{i,j)eE (22) 

This establishes the non-convex DF model. The convex relaxation, CDF, is obtained by relaxing (22) to 
an inequality [6], 

\Sij\ 2 < WiLij V(i,j) eE (23) 

Notice that constraint (23) is also a widely supported second-order cone constraint. 

The complete CDF relaxation of (AC-PF) is presented in Model 3 (CDF-PF). The constraints for 
KCL (12e), generator output limits (12c), and line flow limits (12d) are identical to the AC-PF model 
and constraints (20a) for the voltage bounds are identical to the (SOC-PF) model. Constraints (24c)- 
(24d) capture line power flow in the L-space. Constraints (24a)-(24b) capture the phase angle difference 
constraints, along the lines of (11). Finally constraints (24e) strengthen the relaxation with the power- 
based second-order cone constraint. 


4.3 Equivalence of the Relaxations 

Interestingly, [23] showed that the (SOC-PF) and (CDF-PF) relaxations are equivalent. The key insight 
is that a bijection between the solution sets of (SOC-PF) and (CDF-PF) can be developed as follows. 
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Model 3 The CDF Relaxation of AC Power Flow (CDF-PF) 


variables: 


Sf VieN 
Wi Vi e N 
Lij V(i,j) 6 E 
Stj V(i,j) £ A u E r 

subject to: 

(12c), (12d), (12e), (20a) 


3 (Wi - Z*S ZJ ) V tan(0A)K (W, - Z*^) V(i,j) e A (24a) 

3 (Wi - Z(S,j) > tan(-fl£)R (Wi - Z*S ZJ ) V(i,j) £ A (24b) 

<S l jj + = ZijLij (i,j) £ A (24c) 

Wi - IF, = (Z£% + Z tJ S*) - Z, (i,j) c. E (2 Id) 

\Sn\ 2 ^WiUj V(i,j)eE (24e) 


Given a solution to the (CDF-PF) (Wi, Lij, Sij), assign the (SOC-PF) variables, 

Wij =Wi — ZfjSij ( i,j ) £ E (25a) 

Given a solution to the (SOC-PF) (Wi,Wij, Sij), assign the (CDF-PF) variables, 

= \Yij \ 2 (Wi - Wij — W* + Wj) ( i, j) £ E (26a) 

Note that W and S remain the same in both models. In each of these assignments, all of the model 
constraints are satisfied, demonstrating the bijection. See [23, 3] for a detailed proof. 


5 Extensions of the Power Flow Models 

In the interest of clarity, AC Power Flows, and their relaxations, are most often presented on the purest 
version of the AC power flow equations. However, transmission system test cases include additional 
parameters such as bus shunts fF s ), line charging (fo e ), and transformers (T), which complicate the AC 
power flow equations significantly. This section extends the results of the previous sections to include 
these additional parameters. 

We begin with the definitions of the network’s physical properties. Kirchhoff’s Current Law (KCL) 
is extended to incorporate bus shunts as follows, 

If _ if _ YiVi = ^ lij Vi e N (27) 

(ij)e£uB 8 

Ohm’s Law in extended to include transformers and line charging, 

lij = (Yij + i^j ^ - YjVj V(i, j) £ E (28) 

and the definition of AC power now includes the voltage transformation, 

Sij = ^I*j V(i,j)sE (29) 
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Model 4 The Extended AC Power Flow Feasibility Problem (AC-E-PF) 


variables: 

Variables from (AC-PF) 

subject to: 

(12a),(12b),(12c),(12d),(12g) 


C 9 (. 

|E| 

2 = Yt yieN 

(31a) 



(i,j)eEvE R 


Si 3 = 

(v-4) 

IV I 2 ViV* 

1 l| Y* 3 (i i) e E 

Wil 2 lJ T* 

(31b) 

Sji = 

(v- 4 ) 

\V 3 \ 2 -Y*^ (i,j) e E 

-Llj 

(31c) 


Combining these three properties yields the extended AC Power Flow equations, 

i2 


S? - Sf - (V s 


I VI I 


= E , 

(i,j)eE'uE R 


"V* 

i b 

i |E| 2 

v *ViV; 

*ij 

2 y 

1 l^| 2 

*^3 rp* 

ij 

*ij 

2 y 

) IE I 2 - 

*V*E 

*^3 rp 
-Llj 


■ (i,j)eE 
(i,j) e E 


(30a) 

(30b) 

(30c) 


The complete Extended AC Power Flow Feasibility Problem (AC-E-PF) is presented in Model 4. The 
operational constraints (12a)-(12d) and (12g) remain the same as (AC-PF). Constraints (31a) capture 
the extended KCL and constraints (31b)-(31c) capture the extended Ohm’s Law. 

It is useful to observe that even in this extended formulation, the product of voltages can be factored 
in terms of the V and S variables, as done previously in (10), 

ViV* = Z*T* ( (y* - ^ - Si.) (i, j) e E (32) 

This factorization is useful for extending properties such as the PAD constraints, along the lines of (11), 
and for proving the equivalence of various relaxations. 


5.1 Extensions of the General Properties 

Next the general network properties from Section 3 are redefined in this extended AC power flow context. 


Absolute Square of Ohm’s Law From definition (28), the absolute square of Ohm’s law is, 


resulting in, 

‘ |r «' a (Hf ' 



rri* 

ij 





btMSij) 


(33a) 


(33b) 


A complete derivation appears in Appendix A. 
















Absolute Square of AC Power From definition (29), the absolute square of the definition of ac 
power is, 


c.. q* _ [ ^ r* ) 

Oijdi] — ^ Ji* 1 iJ J 

|S„P - jgfeP 



(34a) 

(34b) 


Absolute Square of Voltage Products The absolute square of voltage products is, 



\ViV*\ 2 = \Vif\Vj\ 2 


(35a) 

(35b) 

(35c) 


The derivation follows similarly to the previous section. Interestingly, the transformer constants cancel 
and this property is unaffected by the model extensions. 


Line Losses From the definitions (30b)-(30c), the power loss on line (i,j) is, 


Sii + Su = V 


m 2 

\Tii \ 2 


rp* 

ij 





(36a) 


Voltage Drop From the definitions (30b)-(30c), the extended voltage drop property is derived by 
subtracting the power on each end of a line (i, j), solving for the |V| 2 , |V)| 2 variables, and eliminating 
the Sji term. 


i C 

Q. . _ Q.. — I V* _ A 


IT l 2 

|r y | 2 


-ITI 2 -V 


WVf 

rp* 

ij 


V? Vi 


resulting in, 

(1 - 3(Zy)&&) S - ITf = ( z foSij + ZijSfi) - ( |Fi|2 


|Iy| a 


Ta I 2 


rp* 

ij 


V*Vi 


\Vj\ 


(37a) 


(37b) 


A complete derivation appears in Appendix A. 


Equivalence of Line Flow Formulations Despite these various extensions, the line loss and 
voltage drop constraints still provide an alternate formulation of the power flow constraints (30b)-(30c) 
namely, 


(30b)-(30c) ^ (36a), (37b) 

which follows from the fact that these two sets of constraints are simply linear combinations of one 
another. 


5.2 The Extended SOC Relaxation 

The Extended SOC relaxation of (AC-E-PF) is presented in Model 5 (SOC-E-PF). The relaxation is 
nearly identical to the simple version as the voltage variables V and the voltage product property are 
unaffected by the model extensions. Simply lifting (AC-E-PF) into the VF-space and strengthening with 
the original second-order cone constraint (20e) completes the relaxation. 
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Model 5 The SOC Relaxation of Extended AC Power Flow (SOC-E-PF) 


variables: 

Variables from (SOC-PF) 

subject to: 

(12c), (12d), (20a), (20b), (20e) 


Q9 < 

if - (*rm 

= Yt yieN 

(38a) 



( i,j)eEuE R 


Sij = 

i 

<s>. 

1 O' 1 

Wi y* Wij / • 771 

(38b) 

Sji = 

(u-4) 

w* 

Wj — Y*- (i,j) g E 

J-ij 

(38c) 


5.3 The Extended Convex DistFlow Relaxation 


Thus far, applications of this DF model have typically focused on distribution systems and, to the best 
of our knowledge, the DF model has not been extended to capture line charging and transformers. Using 
Section 4.2 as a guide, this section derives an extended DF model for AC transmission systems featuring 
bus shunts, line charging, and transformers. 

First, we replace the line flow equations (30b)-(30c) with their alternate formulation based on line 
loss (36a) and voltage drop (37b). Second, we lift the model into the L-space, as follows, 


Wi = \Vi\ 2 i e N 


Lij = |T, | 2 



ViV* 

rji* 

ij 


YIXl 

Tij 






(39a) 

(39b) 


Third, use the absolute square of power property (34b) to strengthen the L-space relaxation, 

\Sa \ 2 = ^2 L ‘i V(i,i)^ (40) 

This establishes the non-convex extended DF model. The convex relaxation, extended CDF, is obtained 
by relaxing (40) to an inequality, 

|Sy| 2 ^b 3 v( i,i)eB (41) 

Notice that constraint (41) is a second-order cone constraint as \Tij\ 2 is a constant. 

The complete extended CDF relaxation of (AC-E-PF) is presented in Model 6 (CDF-E-PF). The 
constraints for generator output limits (12c) and line flow limits (12d) are identical to the AC-PF model 
and the constraints for KCL (38a) and (20a) the voltage bounds are identical to the SOC-E-PF model. 
Constraints (42c)-(42d) capture line power flow in the L-space. Constraints (42a)-(42b) capture the 
phase angle difference constraints, utilizing (32). Finally constraints (42e) strengthen the relaxation 
with a second-order cone constraint based on the extended absolute-square of power property. 


5.4 Equivalence of the Relaxations 

Despite the significant increase in complexity, the (SOC-E-PF) and (CDF-E-PF) are equivalent convex 
relaxations of (AC-E-PF), as demonstrated in this section. 

Theorem 5.1. (CDF-E-PF) is equivalent to (SOC-E-PF). 
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Model 6 The CDF Relaxation of Extended AC Power Flow (CDF-E-PF) 


variables: 

Variables from (CDF-PF) 

subject to: 


(12c), (12d), (20a), (38a) 


tan(-0£)R ( Z*T* ( ( Y* - i-| 


b?A Wi 

T, I 2 K 


< 


3 ( Z*T* ( ( Y* - 


tan(0£)H ( Z*T*i ( ( Y* - 


Sij + Sji — Zjj I Lij + + 

W, 


y u 
Wj a 
Tij I 2 * 

b 


Wj 

IT,;, I 2 ,J 


b C - 


■y 

|Til 2 j 2 


Wi c 

IT;-I 2 ’ J 


I 


V(i,j)e£ 


V(i,i)e£ 




(1 - 3(Zy)&? 


'Til 2 


IE, = (Z*.% + Z,,«v 


l^f 


L ij + KMSij) + 


w 2 

I Til 2 


(i,j)eE 


|5 |2 ^ Wj 
,J I ' IT, I 2 y 


Ty V(i, j) e £ 


(42a) 

(42b) 

(42c) 

(42d) 

(42e) 


Proof. The proof is similar in spirit to those presented in [23, 3] and develops a bijection to establish 
that any solution of one model can be mapped to a solution of the other model. 

We begin by observing that the constraints (12c), (12d), (20a), (38a) are included in both models 
and need not be considered. Hence, this result focuses on the phase angle difference, power flow, and 
second-order cone constraints in each model. Noticing that both models share the variables Wi, Sij, Sji, 
the result assumes these values are unchanged when mapping solutions between the to models. 

(CDF-E-PF) => (SOC-E-PF) 

That is, every solution to (CDF-E-PF) is a solution to (SOC-E-PF). This section will demonstrate that 
given a solution to the (CDF-E-PF) Wi, Lij, Sij, Sji the assignment of the (SOC-E-PF) variables as 
follows, 

Wij - Z*jT*j ((V;* - - Sy) (i,j) e E (43a) 

yields a feasible solution to the (SOC-E-PF) model, i.e. all of the model constraints are satisfied. The 
constraints of interest are the phase angle difference constraints (20b), the power flow constraints (38b)- 
(38c) and the second-order cone constraints (20e). 

The phase angle difference constraints (20b) are satisfied because mapping function (43a) is simply the 
lE-space version of the linear relation (32), which was originally used to transform (20b) into (42a)-(42b). 

It is easy to see that the power flow constraints (38b) are satisfied. Simply substituting (43a) into 
(38b) yields a tautology. The constraints (38c) are more interesting. First we must observe that the line 
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loss and voltage drop constraints, (42c)-(42d), jointly ensure the following property, 


ZijSti = 1^5 - W i + Z * S ^ + ib f 


z o^Tl ~ Z * W ’ I ' 


(44a 


A complete derivation is presented in Appendix B. Combing (43a) and (44a) with the definition of 
constraint (38c) yields a tautology as follows, 


Sa = I 1 Wi - 

J-ij 

Sii W, - S ( ZuTii ( ( V; + i^L 
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_ i LJ V* 7 

L ^ -LinZJi 


Wi 


_ 7 C* 

2 |jr_|2 _r 


+ V-• 
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« • 2 
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b?i 


_ ^ \r * 7 

2 1 2 X< -? ^ 


- Wi + ZTiSa + i- 


Mi 


L *J I 


Wi 

ITijl 2 


Z^ + Ziw, 


(45a) 

(45b) 

(45c) 

(45d) 

(45e) 


Thus demonstrating that the constraints (38c) are satisfied. 

Finally we consider the second-order cone constraints (20e). We begin by developing the value of 
|Wijj 2 , given the variable assignment we have selected as follows, 
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- I 777~ TijZijSij - T V Z* ( i 


bfj W 


2 \Tij | 2 
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|Wj| 2 = {TiiftZiiflSi; 


ZijSijWi — ZijSijWi , 
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(46a 


(ZijWi - Z*Wi + \Tij\ \Zij\ S* - |Tjjn-Zij \ Sij) 


,b c Wi 

1 2 Tj | 2 


it i 2 iz i 2 r c Wi 

1 * j| 1 y| U | Tij | 2 


(46b) 


|Wij | 2 = |Tij| 2 |Zij| 2 |Sij | 2 - (Z*S tj + Z „S*;.W, + 


W 2 

|Tij| 2 


+ (2iSr(Zij)W - 2i\T ij \ 2 \Z ij \ 2 3(S ij )) i 


b c Wi 
2 lTj| 2 


,2,„ |2 fb c W, 


2 | Tj | 2 


(46c) 


|Wj | 2 = Wi |Z, 


2 ( \Tij\ 2 \S, 


Wi 


121 + 3(Sij)b c + {- 


Wj 

|Tij| 2 


- (ZfjSij + ZijSfj) + (1 - 9(Z,j)b c ) 


W 

|Tj| 2 


(46d) 


From this point, we apply the second-order cone constraint from (CDF-E-PF), i.e. (42e), to eliminate 
the |Sij | 2 variable and then apply the voltage drop equation (42d) to simplify the large expression to Wj, 


IWjr -< Wi ( I Zij I' ( Tij + Q(Sij)b c + ( — 


|Wij | 2 s: WiWj 


\Tij \ 2 


Wi ^ — (Z*Sij + ZijSfj) + (1 — Q(Zij)b c ) Wi 


Tij | 2 


(47a) 

(47b) 


Demonstrating that constraint (20e) is satisfied in (SOC-E-PF). With all of the constraints in (SOC-E-PF) sat¬ 
isfied, the proof continues by performing a similar analysis in the reverse direction. 
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(SOC-E-PF) =» (CDF-E-PF) 

That is, every solution to (SOC-E-PF) is a solution to (CDF-E-PF). This section will demonstrate that 
given a solution to the (SOC-E-PF) Wi, Wij, Sij, Sji the assignment of the (CDF-E-PF) variables, 
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\ T a I 2 


Wil 

y* 

ij 


w* 


+ Wi - 


\ 1 

fbfj\ 

) -1 

K 2 J 


Wi 

\Tij\ 2 


V(i,j)eB 


(48a) 


produces a feasible solution to the (CDF-E-PF) model, i.e. all of the model constraints are satisfied. The 
constraints of interest are the phase angle difference constraints (42a)-(42b), the line loss and voltage 
drop constraints (42c)-(42d) and second-order cone constraints (42e). 

The phase angle difference constraints (42a)-(42b) are stratified because the constraint (38b) is simply 
the IF-space version of the linear relation (32), which was originally used to transform (20b) in to (42a)- 
(42b). 

Continuing with line loss constraint (42c), we simply expand the value of Lij as defined in (48a) as 
follows, 


Sji Zij 
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(49a) 

(49b) 


Noticing that this is the extended line loss property (36a) in the W- space demonstrates that the con¬ 
straints hold. Next, the voltage drop constraint (42d) is expanded in a similar way, 


(1 - WZ,,)l,?j) -Wj = (ZfjSij + ZijS *) 

(1 - Q(Zij)btj) -Wj = (.ZfjSij + ZijSfj) 
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■ w, 


(50a) 

(50b) 


Noticing that this is simply the extended voltage drop property (37b) in the TF-space demonstrates that 
these constraints hold. 

Finally we consider the second-order cone constraints (42e). We begin by developing the value of 
|S'ij| 2 , given the model constraints as follows, 
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(51a) 


(51b) 


(51c) 


From this point, factor Wi/ I Tj | 2 from all of the terms. However, to extract this term from \Wij I /I Tij \ 
we apply the SOC constraint from (SOC-E-PF), i.e. (20e), which turns the equation into an inequality. 
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Lastly the assignment of Lij (48a) is used to complete the deriviation as follows, 
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(52a) 

(52b) 

(52c) 


Demonstrating that constraint (42e) is satisfied in (CDF-E-PF), and completing the proof that these two 
extended relaxations are equivalent. □ 


6 Conclusion 

This paper has significantly increased the applicability of the established Convex-DistFlow (CDF) re¬ 
laxation of the power flow equations by proposing an extended CDF model, which is applicable to 
transmission systems. Additionally, it was shown that this extended CDF model defines the same convex 
set as the well known SOC relaxation for transmission systems. Given that these relaxations define the 
same set of solutions, the natural frontier for future work is to conduct a detailed numerical study of 
these two relaxations. Although the relaxation quality may be identical, significant differences in the 
formulations may lead to variations in the computation time of the numerical methods used to solve 
these equations. 
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A Complete Derivation of the Extended Properties 

This appendix provides detailed derivations of a variety of general power flow properties in the extended 
AC-PF model, which are proposed in Section 5.1. Note that these derivations make use of the following 
properties of complex numbers, 

A' + A* = 2K(X) = 2K(X*) 

X - X* = i2Q(X) = -i29f(X*) 
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Extended Voltage Drop 
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B Combination of Line Loss and Voltage Drop 

This appendix provides a detailed derivation of the ZijSfj property, which holds in (CDF-E-PF). This 
property is used in the proof of Section 5.4. In the presence of the line loss equation (42c), the following 
property holds, 
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When combined the voltage drop equation (42d), yields, 
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C The Extended Convex DistFlow in Real Numbers 


In the interest of clean proofs, this document focuses on a complex representation of the extended CDF 
model. In practice however, the a real number representation is often needed for implementation. This 
section presents the real number version of (CDF-E-PF) to aid the implementations of this model. 

First we develop a real number short-hand for Zij Tij as follows, 

F (.Zyly) = t;Zy = ry-ty - *yty (i,j) E E (58a) 

(ZijTij) = tz-j = ryty + ®yty (i,j) e E (58b) 

Notice that the superscripts R, I are used for the real and imaginary parts of a complex number respec- 

tivly. Expanding (32) using this short-hand into real and imaginary components 3l(ViV*),Q(ViV*) we 

have, 

K(W/) = (t* - tz-j |r - tz*pij - tz-jqij (■ i,j ) e E (59a) 

Zs{ViV*) = ^-ty - tzy (%)) ^ + tz ijPH (Aj) e E (59b) 

both of which are used to implement the phase angle difference constraints in any model with v,p,q 
variables. The complete real number implementation of (CDF-E-PF) is presented in Model 7. 
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Model 7 The Extended CDF Relaxation of AC Power Flow (CDF-E-PF) in Real Numbers. 
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